Status of the MODIS Level IB Algorithms and 

Calibration Tables 


X. Xiong and V.V. Salomonson 
Sciences Exploration Directorate 
NASA/GSFC, Greenbelt, MD 20771 

J. Kuyper and L. Tan 
Science Applications International Corp. 
7501 Forbes Blvd., Lanham, MD 20706 

K. Chiang and J. Sun 

Science Systems and Applications, Inc. 

10210 Greenbelt Road, Lanham, MD 20706 

W.L. Bames 

JCET, University of Maryland, Baltimore County 
1000 Hilltop Circle, Baltimore, MD 212250 


Abstract - The Moderate Resolution Imaging Spectroradiometer 
(MODIS) makes observations using 36 spectral bands with 
wavelengths from 0.41 to 14.4pm and nadir spatial resolutions of 
0.25km, 0.5km, and 1km. It is currently operating onboard the 
NASA Earth Observing System (EOS) Terra and Aqua satellites, 
launched in December 1999 and May 2002, respectively. The 
MODIS Level IB (LIB) program converts the sensor’s on-orbit 
responses in digital numbers to radiometrically calibrated and 
geo-located data products for the duration of each mission. Its 
primary data products are top of the atmosphere (TOA) 
reflectance factors for the sensor’s reflective solar bands (RSB) 
and TOA spectral radiances for the thermal emissive bands 
(TEB). The LIB algorithms perform the TEB calibration on a 
scan-by-scan basis using the sensor’s response to the on-board 
blackbody (BB) and other parameters which are stored in Look- 
up Tables (LUTs). The RSB calibration coefficients are processed 
offline and regularly updated through LUTs. In this paper we 
provide a brief description of the MODIS LIB calibration 
algorithms and associated LUTs with emphasis on their recent 
improvements and updates developed for the MODIS collection 5 
processing. We will also discuss sensor on-orbit calibration and 
performance issues that are critical to maintaining LIB data 
product quality, such as changes in the sensor’s response versus 
scan-angle. 

I. Introduction 

Two Moderate Resolution Imaging Spectroradiometer 
(MODIS) instruments, a Proto-flight Model (PFM) and a Flight 
Model 1 (FM1), were designed and built by Raytheon / Santa 
Barbara Remote Sensing (SBRS), Goleta, California, for the 
NASA’s Earth Observing System (EOS) missions that promote 
the continuity of essential science and environmental data 
records [1-5]. The MODIS PFM onboard the Terra spacecraft 
has been in operation for more than 5 years since its launch in 
December 1999 and the FM1 for nearly 3 years onboard the 
Aqua spacecraft launched in May 2002. Each instrument makes 
continuous observations over a wide field-of-view (FOV) in 36 


spectral bands with wavelengths from 0.41 to 14.4pm and at 
nadir spatial resolutions of 0.25km, 0.5km, and 1km. Table 1 
lists key sensor specified parameters and characteristics. 


Table 1 MODIS Specified Parameters and Characteristics 


Orbit 

705km (altitude) 

Sun-synchronous, near-polar, circular 
10:30AM descending node (Terra) 

1 :30PM ascending node (Aqua) 

Spectral Range 

20 reflective solar bands (1-19 and 26) from 
0.41 to 2.2p, 

16 thermal emissive bands (20-25 and 27-36) 
from 3.7 to 14.4p 

Spatial Resolution 
(nadir) 

0.25km (bands 1-2, 40 detectors along-track and 
4 sub-samples each cross-track) 

0.5km (bands 3-7, 20 detectors along-track and 
2 sub-samples each cross-track) 

1km (bands 8-36, 10 detectors along-track) 

Swath 

10 (1km resolution) pixels along-track by 1354 
(1km resolution) pixels cross-track 

Scan Rate 

20.3 rpm ( 1 .477s / scan) 

Scan Range 

±55° from nadir or angle of incident (AO I) to 
the scam mirror from 10.5° to 65.5° 

Digitization 

12 bits ADC 

Data Rate 

1 1Mbps daytime 
3.3Mbps nighttime 

On-board 
Calibrators (OBCs) 

Solar diffuser (SD) 

Solar diffuser stability monitor (SDSM) 
Blackbody (BB ) 

Spectro-radiometric calibration assembly 
(SRCA) 

Focal Plane 
Assemblies (FPAs) 

Visible (VIS) 

Near infrared (NIR) 

Short- and mid-wave infrared (SMIR) 
Long- wave infrared (LWIR) 

Spectral Bands and 
Detectors 

36 bands 

490 detectors (including time-delay and 
integration detectors used for the high/low gain 
bands 13 and 14) 


The MODIS Level IB (LIB) algorithms convert the 
sensor’s responses in digital numbers (DN) to top of the 
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atmosphere (TOA) reflectance factors for the reflective solar 
bands (RSB) and TOA spectral radiances for the thermal 
emissive bands (TEB). They are designed based on instrument 
performance determined from comprehensive pre-launch and 
on-orbit calibration and characterization [6-9]. Two separate 
versions of the LIB software and corresponding Look-up 
Tables (LUTs), one for Terra MOD1S and the other for Aqua 
MOD1S, are maintained to accommodate differences between 
the two sensors and to track their on-orbit changes 
independently. The NASA MODIS Characterization Support 
Team (MCST) is responsible for the tasks of performing sensor 
calibration and characterization, maintaining and improving the 
LIB algorithms, and deriving and updating the calibration 
parameters in the LIB LUTs. In this paper we provide a short 
description of the MODIS LIB calibration algorithms and 
associated LUTs and highlight the latest improvements and 
updates implemented for the MODIS collection 5 processing. 

II. Overview of MODIS LIB Calibration Algorithms 

The MODIS LIB calibration is performed for each band, 
detector, sub-sample, and mirror side (BDSM). The algorithm 
consists of two major modules: one for the reflective solar 
bands (RSB) and the other for the thermal emissive bands 
(TEB). Details of the MODIS TEB and RSB calibration 
algorithms are provided in papers given by Xiong, et al [10- 
11]. The LIB input includes the L1A data (raw counts), geo- 
location data, and LIB LUT files for the TEB, RSB and quality 
assurance (QA). A simplified LIB data processing flow 
diagram is shown in Figure 1. 


START 


1 READ LOOKUP TABLES I 

— 




2 

< ■ 

o 

to 


PREPROCESS OBC DATA 1 

READ LI A EV SCAN 

1 



OPEN FILES; INITIALIZE DATA 

| EMISSIVE CALIBRATION | 



<^DAY" 

^\NO 


UJ 


A 1 

CALCULATE GRANULE METADATA 

O 

8 

O 

Y ^T 1 

T 

WRITE GRANULE METADATA 

| REFLECTIVE CALIBRATION | 

1 

| AGGREGATE 250M, 500M | 

1 

CLOSE L1A GRANULE 

WRrTE LIB EV SCAN | 

1 



CLOSE LIB GRANULE 



Figure 1 MODIS Level IB (LIB) Flow Diagram 


The thermal emissive bands (20-25 and 27-36) are 
calibrated by an on-board blackbody on a scan-by-scan basis 
using a quadratic algorithm that corrects for instrument 
background (from the sensor’s space view), thermal emissions 
from the scan mirror and scan cavity, and response versus scan 
angle (RVS). An optical leak correction is applied to bands 32- 
36 (Terra only). The retrieved Earth view (EV) radiance is 
converted to a scaled integer (SI). 


The reflective solar bands (1-19 and 26) are calibrated 
using a solar diffuser (SD) and solar diffuser stability monitor 
(SDSM) system. The SD/SDSM calibration is typically 
performed bi-weekly (weekly for the first year) and analyzed 
off line by the MCST. The RSB calibration uses a linear 
algorithm with calibration coefficients updated through the 
LUTs. The algorithm corrects for instrument background, 
temperature, and response versus scan angle (RVS). A thermal 
leak correction is applied to the SWIR bands (5-7 and 26). The 
retrieved EV reflectance factor is converted to a scaled integer 
(SI). The RSB calibration is executed after the TEB calibration. 

The LIB output consists of calibrated radiances and 
reflectance factors and is organized in three Hierarchical Data 
Format (HDF) files corresponding to MODIS’ three spatial 
resolutions and their associated metadata files. The following 
expressions are used to convert the SI to the scene TOA 
radiances or reflectance factors. 

For TEB radiances: 

Radiance = radiance _scales * (SI — radiance _offsets) 

For RSB reflectance factors: 

Reflectance = reflectance _scales* (SI -reflectance ^offsets) 

For RSB radiances: 

_ ,. _ Solar Irradiance over pi 

Radiance = Reflectance = = =f— 

( Earth - Sun _ Distance ) 

The band dependent coefficients of radiance_scales, 
radiance _offsets, reflectance _scales, and reflectance ojfsets 
are stored as scientific data set (SDS) attributes within the 
Level IB output files. The “Earth-Sun Distance” and 

“ Solar _lrradiance_over_pr are provided as global attributes. 
The LIB also produces a separate file containing on-board 
calibration (OBC) data and key instrument telemetry. The 
OBC file is primarily used for off-line calibration and 
characterization purposes. 

III. MODIS LIB Look-Up Tables 

The MODIS Level IB (LIB) code Look-up Tables (LUTs) 
are organized into three groups: one for the TEB calibration, 
another for the RSB calibration, and a third one for quality 
assurance (QA). Each individual LUT is implemented in the 
appropriate HDF file as either a global attribute or a Scientific 
Data Set (SDS) [12], Currently there are 98 LUTs (40 TEB, 34 
RSB, and 24 QA) for the Terra MODIS LIB code and 101 
LUTs (43 TEB, 34 RSB, and 24 QA) for the Aqua MODIS 
LIB code. Aqua MODIS has 3 extra TEB LUTs compared 
with Terra MODIS. These extra LUTs are used to provide pre- 
determined linear gain coefficients for bands 33, 35, and 36 at 
high BB temperatures during BB warn-up process when their 
responses saturate (a feature that only appears in the Aqua 
MODIS). 



Some of the LUTs included in the HDF files are not used 
for calibration, such as those used for serial numbers, metadata 
items, and code version numbers. Some of the LUTs are 
constant and some may depend on the data collection time in 
either a “step function’’ or “piecewise linear” style. Among 
calibration LUTs, some were derived from pre-launch 
calibration and characterization, such as the relative spectral 
response (RSR), the temperature correction coefficients, and 
corresponding reference temperatures. Some calibration LUTs 
are derived from on-orbit measurements and updated as 
needed. Examples include the RSB response versus scan angle 
(RVS) and SD calibration coefficients for the reflectance factor 
retrieval. The MCST at NASA/GSFC is responsible for 
determining the contents of all LIB LUTs and for maintaining 
and updating (adding, deleting, or changing) them. 

IV. Status OF MODIS LIB and LUTs 

Terra MODIS has been in operation for more than five 
years since its launch in December 1999. Thus far there have 
been four major configurations used for the instrument 
operation. As shown in Table 2 that there are 12 versions of the 
LIB code used in data production. Each change was made to 
improve the LIB algorithms or to mitigate code deficiencies. 
Many versions of the LUTs have been updated either to the 
track sensor’s response changes over time or to accommodate 
changes in operational configuration. 


Table 2 Terra MODIS LIB code and LUT versions 


Code Version 

Production Time 

LUTs Version 

V2.3.2 

03/17/2000 (077 2000) 00:00 

3 

V2.4.2 

06/19/2000(171 2000) 00:00 

0 

V2.4.3 

08/18/2000(231 2000) 14:00 

1 

V2.4.4 

10/13/2000 (287 2000) 19:55 

0-2 

V2.5.4 

11/23/2000 (328 2000) 15:55 

0 

V2.5.5 

02/13/2001 (044 2001) 13:55 

1,2 

V3.0.0 

05/24/2001 (144 2001) 00:00 

2-7 

V3.0.1 

02/25/2002 (056 2002) 00:00 

0, 1 

V4.1.2 

01/30/2003 (030 2003) 01:55 

1-11 

V4.2.0 

08/22/2003 (234 2003) 02:00 

3-9 

V4.3.0 

12/22/2003 (356 2003) 22:35 

1-9, 11-25 

V5.0.6 

03/07/2005 (066 2005) 23:55 

1,2 


Since its launch in May 2002, Aqua MODIS has been in 
operation for nearly three years. Unlike its Terra predecessor, 
the Aqua MODIS has been operated using the same B-side 
electronic configuration for the entire mission. Table 3 is a list 
of its LIB code and LUT versions. 


Table 3 Aqua MODIS LIB code and LUT versions 


Code Version 

Production Time 

LUTs Version 

V3.1.0 

06/07/2002 (158 2002) 18:10 

0-3 

V4.1.1 

10/31/2002 (304 2002) 00:15 

0,1 

V4.1.3 

01/22/2003 (022 2003) 09:55 

0-11 

V4.2.1 

08/21/2003 (233 2003) 12:00 

4-8 

V4.3.1 

06/02/2002 (183 2002) 15:25 

1-16 

V5.0.7 

Scheduled for May 2005 

1 


The latest MODIS data collection will be produced using 
version 5 of the LIB code. The Terra collection 5 (V5.0.6) was 
started in March 2005 and the Aqua collection 5 will start in 
May 2005. In addition to the NASA Goddard Distributed 
Active Archive Center (GDAAC), selected LIB products have 
also been produced at the MODIS Adaptive Processing System 
(MODAPS) and at the MODIS Ocean Color Data Processing 
System. 

Only a few minor code changes were made in LIB version 
5, including improving code portability, enabling B21 
calibration with mirror side dependency (for reducing image 
striping), and adding new LUTs for SWIR correction with 
flexible detector dependency in order to avoid using noisy 
detectors in the correction algorithm. Overall improvements are 
primarily coming from better calibration and LUTs that are 
developed using long-term consistent calibration data. 

In version 5 the Terra TEB will use the response versus 
scan-angle (RVS) derived from on-orbit spacecraft pitch 
maneuvers made on March 26 and April 14, 2003 for the entire 
mission, replacing the relative RVS determined from scanning 
the closed nadir aperture door. For the reflective solar bands of 
both Terra and Aqua MODIS, on-orbit RVS changes of mirror 
side one are updated using each sensor’s response trending 
from SD calibrations and lunar observations. The RSB RVS 
changes of mirror side two are derived relative to mirror side 
one using the trending of mirror side response ratios at various 
angles of incidence (AOI). Several QA LUTs have been 
updated based on the detectors’ on-orbit performance. 
Complete histories of LIB code versions and LUT updates for 
both Terra and Aqua MODIS are posted on the MCST web 
page ( http://www.mcst.ssai.biz/mcstweb/LlB/ product. html ). 

Other documents that describe the algorithms, function and 
applications of the LIB code and LUTs include “MODIS Level 
IB Algorithm Theoretical Basis Document”, “MODIS Level 
IB Product User’s Guide”, “MODIS LUT Information Guide”, 
“MODIS Level IB In-Granule Calibration Code High Level 
Design”, and “MODIS Level IB Products Data Dictionary”. 
These are also provided at the MCST web page (home page 
http://www.mcst.ssai.biz/mcstweb/ ). 

V. Summary 

Since their launches in December 1999 and May 2002, 
Terra and Aqua MODIS together have produced nearly eight 
years of data products for the science community’s studies of 
the Earth system and its land, oceans, and atmosphere. 
Working together with the MODIS Science Team members 
and sensor vendor representatives, the NASA MODIS 
Characterization Support Team (MCST) is responsible for 
providing instrument on-orbit calibration and characterization, 
maintaining and improving the LIB algorithms, and updating 
the calibration parameters in the LIB LUTs. Currently two 
separate versions of the LIB software and LUTs have been 
maintained and updated for the data processing, one for the 
Terra and the other for the Aqua MODIS. This paper has 
presented the status of the MODIS LIB algorithms and 
associated LUTs, including the latest improvements and 
updates made for the MODIS collection 5 processing. 
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